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Abstract: We present a construction kit for calculating two- loop beta functions 
in N = 1 supersymmetric theories for the operators of the superpotential using 
supergraph techniques. In particular, it allows to compute the beta functions for 
every desired, even higher dimensional, operator of the superpotential from the 
wavefunction renormalization constants of the theory. We apply this method to 
calculate the two-loop beta functions for the lowest-dimensional effective neutrino 
mass operator in the Minimal Supersymmetric Standard Model (MSSM) and for the 
Yukawa couplings in the MSSM extended by singlet superfields and the mass matrix 
for the latter. Our method can be applied to any N = l supersymmetric theory. 
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1. Introduction 



In order to compare experimental results with predictions from models beyond the 
Standard Model (SM), like unified theories, it is essential to evolve the parameters of 
the models from high to low energies. This is accomplished with the renormalization 
group equations (RGE's) for the operators in the theory. Besides the renormaliz- 
able operators there can be higher dimensional, non-renormalizable operators if the 
theory is considered as effective. 

In a previous study [1], we derived a general method for calculating /3-functions 
from counterterms in MS-like renormalization schemes, which works for tensorial 
quantities. It is simplified considerably in supersymmetric (SUSY) theories, since 
due to the non-renormalization theorem [2,3] only wavefunction renormalization has 
to be considered for operators of the superpotential. However, in a component field 
description, no use can be made of the theorem with respect to gauge loop corrections 
since it is no longer manifest when a supergauge, as for example Wess-Zumino-gauge, 
has been fixed. The supergraph technique [4-7], on the other hand, allows to use 
the non-renormalization theorem, since SUSY is kept manifest. Moreover, it has the 
advantage that the number of independent diagrams is clearly reduced compared to 
the component field calculations. 

We therefore present a method to calculate /^-functions in supersymmetric the- 
ories for operators of the superpotential from wavefunction renormalization. These 
operators may be non-renormalizable since for the latter the non-renormalization 
theorem holds as well [8], and they do not affect the wavefunction renormalization 
constants in leading order in an effective field theory expansion. As an application, 
we consider the Minimal Supersymmetric Standard Model (MSSM) extended by 
singlet superfields, which contain right-handed neutrinos relevant for models of neu- 
trino mass. We compute and specify the wavefunction renormalization constants. 
From these, the two loop RGE's for the Yukawa couplings and for a possible mass 
matrix for the singlet neutrino superfields are obtained by the supergraph method. 
The technique can be used to compute the two loop RGE's for every desired higher 
dimensional operator of the superpotential, since the non-renormalization theorem 
guarantees that no vertex corrections contribute. Furthermore, we consider the low- 
est dimensional neutrino mass operator. Its RGE is known at the one loop level 
for the SM [9,10,1], Two Higgs Doublet Models [10,11] and for the MSSM [9-11]. 
With the supergraph method we calculate the two loop RGE for the neutrino mass 
operator from the MSSM wavefunction renormalization constants. 
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2. Review of the Supergraph Method 



Consider a general supersymmetric gauge theory with a gauge part described by the 
usual Lagrangian 

n,,,, r. „ s 



% = J d 4 * I ex P( 2 S • * 0) + \ J d2 ° XJ W «W na + h.c. 

~l~ =^Ghost -S^Gauge Fixing 



d 2 0)T + h.c. 



n=l 



where 



- 1 - D 2 [exp(2 5 „ V") £> Q exp(-2 5n V n )] 

8^n 



5" dim Gn 

5-V:=^5nV" and V" = ]T 
n=l A=l 

The renormalizable part of the superpotential reads 



' n n 



6 E A fo - fc) #^*W*«. 



(2.1) 

(2.2a) 
(2.2b) 

(2.3) 



i,j,k=l 



Possible mass terms are ignored for the present as they do not affect the /3-functions 
of the model. The N® superfields <fl>W transform under the irreducible representa- 
tions (irreps) R$ ) x • • • x R,f of the gauge group G x ® ■ ■ ■ ® G s . {(T^)]}^™ Gn 
denote the generators of G n . The indices . . . run over all irreps, families and the 
representation space. 

2.1 Wavefunction Renormalization Constants 

Due to the non-renormalization theorem, the RGE's for operators of the superpoten- 
tial are governed by the wavefunction renormalization constants for the superfields 

(i) 

Zij = Ijj + 5Zij which relate the bare <E>g and the renormalized superfields, 



A'* 



(2.4) 



Results are obtained with dimensional regularization via dimensional reduction [12,13]. 
At the one loop level and in d = 4 — e dimensions, 5 Zij is given by 

N® S 
/2 Ku X M "4^fln C 2 (ll®) Sij 



sz 



(1) 



1 1 

(4^7 



(2.5) 



k,£=l n=l 

In equation (|2.5|) and in the following, we use the group-theoretical constants 

C1 (G)5 AB :=Y,f ACD f B CD, (2.6a) 



C,D 



C2(R) Sab 
i{R) 5 AB 



Y,(T A V 



Ti(J A T B ) 



(2.6b) 
(2.6c) 
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with the matrix representations {T^}^ 1 ™^ of the generators of G corresponding to 
the irrep R and the structure constants f bc- £(R) is known as Dynkin index of 
the irrep R and C2(R) as the quadratic Casimir. They are related by 



c 2 (R) 



dimG 
dimi? 



e(R), 



(2.7) 



with dimCr and dimi? being the dimension of the group G and the irrep R, re- 
spectively. Often the generators of the irrep N of SU(iV) are normalized such that 



(AT") = \ holds. C2 can then be obtained via C2(N) 



N*-l 



while for a U(l) theory 



both £(R) and c 2 (R) are replaced by q 2 where q is the U(l) charge of <D>. For any 
non-trivial irrep of SU(iV) the invariant ci(N) is given by N. 

On the two-loop level the renormalization group equations are determined from 
the formula [14] 



-SZ, 



(2) 



-2 + e 
(4vr) 4 e 2 



4 £ g 2 n c 2 (R$) g 2 m c 2 (R$) b %3 

n,m=l 
S 

+ 2 ^9^2 (4°) & " 3 Ci (G n )) % 

71=1 

S AT* 

n=l fc,£=l 
AT* 

~~ « X/ Kkl^tst Kst^jkr 



fe,i,r,s,t=l 



(2i 



Here £ n is defined by 



/ •- 



i= i dim [R%' 



(2.9) 



In addition to the superpotential of equation (|2.3| ) , higher dimensional operators 
may appear in the superpotential of an effective theory. These operators are gener- 
ally suppressed by inverse powers of a large mass scale Mx- Though these operators 
are non-renormalizable by power counting, in the effective field theory approach one 
can renormalize the theory in an expansion in inverse powers of Mx- In the leading 
order of this expansion, the higher dimensional operators do not contribute to the 
wavefunction renormalization. 

In supergraphs, we represent chiral superfields as straight double lines while 
vector-superfields are indicated by wiggly double lines, 
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The diagrams relevant for the calculation of the wavefunction renormalization con- 
stants for the matter superfield, equation (|2.5|) and (|2.8| ), are shown in figure [j] and 
||, respectively. 




(g) (h) (i) 

Figure 2: Two-loop supergraphs which contribute to the cjxft propagator. A blob denotes the relevant 
one-particle irreducible graph including any one-loop counterterm that may be required [14]. 
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2.2 /^-Functions from Wavefunction Renormalization Constants 

To calculate /3-functions from the wavefunction renormalization constants, it is con- 
venient to subdivide the general indices into indices {r, s, . . . } for the 
irreducible representations, {/, g, . . . } for the particle families and {a, b, . . . } for the 
representation space, i.e. i = (r,f,a). The wavefunction renormalization constants 
Z{j are diagonal with respect to the representation and the representation space 
indices and are matrices in flavour-space. In the rest of the paper we will write 

= Z r and suppress flavour and representation space indices. 

Due to the non-renormalization theorem, in a supersymmetric theory a bare 
quantity Qb of the superpotential and the corresponding renormalized one, Q, are 
related by 



qb= (n z "' ) Q» DQt (n 

\r€l J \seJ 



z: 



(2.10) 



where matrix multiplication with respect to flavour indices is implicit and the sets 
of superfield indices I = {1, . . . M} and J = {M + 1, . . . N} denote the wavefunction 
renormalization constants multiplied from the left and the right respectively. Q may 
correspond to a renormalizable or non-renormalizable operator. 

The wavefunction renormalization constants Z r can be expanded in e, 



(2.11) 



k>l 



Following the steps of the derivation in [1] , we find 

'dZ« i 



. A 



rel 



dV A 



V a 



dZ, 



r.l 



dV A 



V A 



■Q 



(2.12) 



where {Va} denotes the set of all variables of the theory including the variable Q 
under consideration and f3n is the usual /3-function, defined by 



(2.13) 



Note that in equation (2.12), for complex quantities Va we have to treat the complex 
conjugates F| as independent variables. We use the notation [1] 



dF 
dx 



( dF 

^dF 
^dx~ Vn 

n aXn 

Y—v 



for scalar s x, y 

for vectors x = (x m ),y = (y m ) 

for matrices x = (x mn ),y = (y r , 
etc. 



(2.14) 
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and Dy A is related to the mass dimension of Va as indicated in equation ( 2.1CQ . 
Equation ( [2.12| ) allows to compute the /3-functions directly from the wavefunction 
renormalization constants, calculated with the supergraph technique 1 . It has a form 
that can easily be used for computer algebra calculations. 



3. Applications 

3.1 Two-Loop /3-Functions in the MSSM Extended by Singlet Superfields 

We consider a supersymmetric model containing the same fields as the MSSM and 
additionally the singlet "neutrino" superfield which we will denote by v. In order to 
obtain the /3-functions for the Yukawa matrices and a mass matrix for the neutrino 
superfield, we will omit the soft SUSY breaking terms 2 , since they do not affect 
the considered /3-functions above the scale of the soft supersymmetry breaking mass 
terms. Threshold effects at low energy scales are e.g. discussed in [19]. 
Thus the Yukawa part of the superpotential is given by 

^Yukawa = (Y e ) gf e C ^ 8^ + (Y v ) gf^h® E a H{ 

+(Y d ) gf sf^e a \{ + {Y u ) gf ^\e T f\{ . (3.1) 

The superfields e c , <d c and u c contain the SU(2)L-singlet charged leptons, down-type 
quarks and up-type quarks, respectively, and q contains the SU(2)l quark doublets. 
Their quantum numbers are specified in table |]. In addition we consider a mass 
term for the singlet neutrino superfield 

^Mass = \v°fM fg v c v , (3.2) 

that may be relevant for generating neutrino masses in the see-saw scenario. 

The quantum numbers of the superfields are listed in table |]. Note that we use 
GUT charge normalization for the U(l)y charge. 



Field 


lh« 


h( 2 ) 


q 






e c 






i 

2 


+\ 




+ § 


2 
3 


~\ +1 





SU(2) L 


2 


2 


2 


1 


1 


2 1 


1 


SU(3) C 


1 


1 


3 


3 


3 


1 1 


1 



Table 1: Quantum numbers of the superfields. qY denotes the U(1)y charge in GUT normalization. 



Another way of calculating superpotential /3-functions is based on the superfield anomalous 
dimensions, which have been calculated to the 3-loop level in [15]. 

2 The calculation of the RGE's for the soft SUSY breaking operators can be found in [16-18]. 
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Using equation fl2.5|) , this leads to the 1/e-coefficients of the wavefunction re- 
normalization constants 

-(4n) 2 Z$ )A = 6 Tr (yj ■ Y d ) + 2 Tr (yj ■ Y e ) - -gf - 3g 2 , (3.3a) 

-(4n) 2 = 6 Tr (yj . Yu ) + 2 Tr (yj . y v ) - ~ g\ - 3g 2 , (3.3b) 

-(4vr) 2 Z§ = 2 Fj • Y d + 2Y% ■ Y u - 1 g\ - 3g 2 2 - y <? 2 , (3.3c) 

-(4«) 2 Z$ tl =4Y d *.Y d T -±g 2 - 1 l g 2 , (3.3d) 

-(4vr) 2 ZW, = AY U * ■ Yt ~ f 5 9l ~ Y 9l , (3-3e) 

-(4vr) 2 Z$ = 2 yj ■ Y e + 2 yj ■ Y v - ^g 2 - 3 5 | , (3.3f) 

-(4^) 2 z« 1 = 4y;.yJ-^?, (3. 3g ) 
-(4^^ = 4y;.yJ, (3.3h) 

where the Y's as well as the last six Z-factors are of course matrices in flavour space. 
From the two-loop diagrams we get 

Z (2) 

-(4vr) 4 = -9 Tr(Y d ■ yj ■ Y d ■ yj) - 3 Tr(Y u ■ Yj ■ Y d • yj) 

-3 Tr(y e • yj • y . yj) - Tr(y • yj . y„ . yj) 

~ g j TV(yj • y d ) + jj 5 2 Tr(y e t • y e ) + 16 5 2 Tr(yj • y d ) 

207 4 9 „ 9 15 4 . , . 

+ Too^ + Io 5 ^ 2 + T 92 ' (3 ' 4a) 

Z (2) 

-(4vr) 4 = -3TV(y u • yj • Yd ■ yj) - 9Tr(y u • yj . y u . yj) 
-TY(y, • yj • y • yj) - 3 tf(y„ • yj • y„ • yj) 

+ 4 - 9 2 Tr(yj • Y u ) + 16 5 2 Tr(yj . Yu ) 
207 4 9 2 2 15 4 ,„ N 

+ ioo^ + Io 5?52 + T 52 ' (3 - 4b) 

Z (2) 

-(4vr) 4 ^i = -2yJ.y ci .yJ.y d -2yj.y 1 .yj.y 1 

-3 yj • Y d Tr(y d • yj) - 3 yj • y u Tt(y u ■ yj) 
-yj • Y d Tr(y • yj) - yj • y u Tr(y, • yj) 

+jjsi^ t -^+ 4 - 5 2 yj-y„ 

199 4 1 2 2 15 4 

H 9i H <?T 99 H 99 

900 yi 10 yia2 4 y2 

+^ff?53 2 +85 2 2 9i-^3, (3.4c) 



- 7 - 



_( 47r )4^i = _ 2Y * . Y J -Y d * -Yj -2Y d * -Yj -Y* -Yj 

-6 Y d * ■ Yj Tr(Y d • yj) - 2 Y d * • Yj Tr(y • yj) 

202 A 32 9984 , , , . 

+ 225 * + 45 *"9*' (3 - 4d) 

_ ( 4^^i = -2y u *.yJ.y;.y u T -2y;.y u T -y ? :.y u T 

-2 y tt * • yj TV(y • yj) - 6 y; • yj Tr(y M • yj) 

-^?^:->; T + 6<? 2 2 y;-y u T 

856 a 128 9984 ... . . 

+ 2l^ + 45-^-9^' ^ 

Z {2) 

-(4vr) 4 = -2y e t • y • y e t . y e - 2yj . y„ • yj • y 

-3 yj • y Tr(Y d ■ yj) - yj • y Tr(y • yj) 
-3 yj • y Tr(y u • yj) - yj • y Tr(y • yj) 

+ ^? y e t -n + ^^ + ^^ 52 2 + 1 | 52 4 , (3.4f) 

-^^ = -2y:.y-.y:-y--2y:.y?. v .y- 

-6 y; • yJ Tr(y • yj) - 2 y; • y e T Tr(y • yj) 

~\ 9l Y: ■ Y e T + &gl Y: • *f + f^r , (3.4g) 

Z {2) 

-(4tt) 4 = -2 y; • Y e T ■ y; • yj - 2 y; • yj • y; • yj 

-6 y; • yj Tr(y u • yj) - 2 y; • yJ Tr(y • yj) 

+\g\ y; • yj + e 92 2 y; • yj , (3.4h) 
5 

respectively. From these, the two-loop Yukawa RGE's are derived, 

^ = 7±~f® + tAz4? , (3-5) 



where x € {d, u, e, z/|. Using equation ( 2.12p , the one-loop contributions to the 
/5-functions are given by 



3 {1) 



Y d ■ I 3 yj • Y d + yj • y u + 3 Tr(yj • Y d ) + TY(yj • y) 

- ttSi - 3#f - ^ 53 1 , (3.6a) 



15 J1 " z 3 
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(3 { y>=Y u .{ yj -Y d + 3Yj -Y u + Tr(Yj -Y v ) + 3TY(yj -Y u ) 



$J = Y e ■ 3 yj ■ Y e + yj • y„ + 3 Tr(yj • y d ) + Tr(y e t ■ y e ) 



= y v ■ yj ■ y e + 3 Yj ■ Y v + 3 TY(yj • Y u ) + Tr(yj • Y v ) 



$l= Y d-\ -^Y}-Y d .Yj.Y d -2Y^Y u .Y}.Y d -2Y^Y u -Y^Y u 



- 9Tr(y d • yj • y d • yj) - 3Tr(y d • yj • y u . yj) 

- 3 Tr(y e • yj • y e ■ yj) - Tr(y e • yj • y v ■ yj) 

- gyj • y,Tr(y d ■ yj) - syj ■ y,Tr(y e • yj) 

- yj • Y u Tr(Y v • yj) - 3 yj • Y u Tr(Y u ■ yj) 

+ 6 g \ yj ■ y d + yj ■ y d + ±gj yj • y u 

-\g\ Tr (yj • y d ) + 1 g\ Tr (yj • y e ) + ie g \ Tr(yj ■ Y d ) 

+ ^gt + 9lgl + ^gi + lahl + zglgi-^gi } , (3.7a) 

/£!? = y « • {" 2y i • y ^ • y i -Y d -2Yj.Y d . yj • y u - 4yj ■ y u • yj • y u 

- syj • y d Tr(y, • yj) - yj ■ y d Tr(y e ■ yj) 

- 9 yj • y u Tr(y u • yj) - 3 yj • y u Tr(y, • yj) 

- 3 TV(y u • yj • y d • yj) - 9 Tr(y u • yj • y u • yj) 

- Tr(Y v • yj • Y e ■ yj) - 3 Tr(Y„ ■ Yj ■ Y v ■ yj) 
+ 2 -g\ yj • Y d + ^ yj . Yu + 6 5 2 y j . y u 

+ 1 5? Tr(yj • y u ) + i6si Tr(yj • y„) + H 

+ <7? 5 2 2 + y 5 2 4 + ^ g\ sH + 8 si g\ - y 53 } . (3-7b) 




(3.6b) 




(3.6c) 




(3.6d) 



and the two-loop contributions are 
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5 • I - 4F e t ■ Y e ■ yj • Y e - 2Yj • Y v ■ yj -Y e -2Yj-Y u - yj ■ Y v 

- 9 yj • y Tr(y • yj) - 3 yj ■ y Tr(y • y e t) 

- yj • y Tr (y • yj) - 3 yj ■ y Tr (y u ■ yj) 

- 9Tr(y • yj • y • yj) - 3Tr(y • yj • y u ■ yj) 

- 3 Tr(y • y e t • y ■ y e t) - Tr(y • yj • y • yj) + \& Tr(yj ■ y ) 
+ 6 g \ y e t • y - 1 5l 2 Tr (yj • y ) + 16 g \ Tr (yj • y ) 

27 4 9 2 o 15 4 1 < . 

+ y5i + 5^52 2 + y52 4 ] , (3.7c) 

y • | - 2Y} ■ y • y e t • y - 2yj • y • yj • y - 4yj • y„ . yj . y 

- 3 yj • y Tr(y • yj") - yj ■ y Tr(y • y e t) 

- 3 yj • y Tr(y • yj) - 9 yj • y Tr(y u • yj) 

- Tr(y • yj • y • yj) - 3 Tr(yj • y • yj • y ) 

- 3Tr(y M • yj • y • yj) - 9Tr(y u • yj • y u • yj) 

+ I 5l 2 yj . y + 6 2 yt . y, + 65 2yt . ^ 
5 5 

+ ~ <?? Tr (yj • y M ) + 16 g\ Tr (yj • y ) 

207 4 9 9 9 15 4 1 , „ 

+ ^o 9l + 5 9l92 + Y 9 ' M • (3 ' 7d) 

Note that the two-loop MSSM RGE's for Y^, Y u and y are easily obtained by setting 
y, = 0. The effort is clearly reduced compared to component field calculations [20- 
22]. 

3.2 Two-Loop /3-Function for the Effective Neutrino Mass Operator 

We now apply our method to calculate the /3-function for the lowest dimensional 
effective neutrino mass operator, which is contained in the .F-term of the superpo- 
tential 

It can e.g. be obtained by integrating out the singlet superfield v of the model 
described in section |3.1| at leading order in the effective theory. The /3-function can 
easily be computed using our method. Substituting D gi = Dy x = \ with i G {1, 2, 3} 
and x G {u, d, e}, we get from equation ( 2, 12] ) 



Pk = - Z h^),\ ■ K ~ \ Z l,l ■ K ~ \ K ■ • ( 3 - 9 ) 
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We can thus write the /3-function for n in the form 

f3 K = X T ■ k + k ■ X + a k 



(3.10) 



where the complete flavour diagonal part is contained in a. We further split X = 
X^ + X^ and a = + into their one loop and two loop part. Plugging in 
the wavefunction renormalization constants of equation ( p.3b| ) and ( |3.3iD and setting 
Y v = 0, our method reproduces the one loop results of [9-11] 



(4vr) 2 X« = y e t • Y e , 
(4vr) 2 «( 1 ) 



--gf-6g 2 2 +6Tr(Yj.Y u ) 



(3.11a) 
(3.11b) 



Note that for U(1)y, we use GUT charge normalization as specified in table |T|. At 
two- loop, with the wavefunction renormalization constants given in equations ( p. 4b ) 
and ( 3.4f ), we obtain 

(4vr) 4 X( 2 ) 



and 



(4vr) 4 a 



2 Y^ • Y e ' Y^ • Y e 
+ (~9i- Tr(F e • y e t) - 3 Tr(Y d • Y*)) Y^ ■ Y e 

-6 Tr(Y u ■ Y\ ■ Y d ■ Fj) - 18 Tr(y M • Fj • Y u ■ yj) 

+ § 0? iv(y u t • y„) + 32 <? 3 2 Tr(yJ • Y u ) 
207 / 18 n 9 j 

+ ^5" 5 ' + y 5 ' ff2+ 52 • 



(3.12) 



(3.13) 



3.3 Two-Loop /3-Function for the Mass of the Singlet Superfield 

From the wavefunction renormalization constants of the MSSM extended by singlet 
superfields given in section |3.1| , the /3-function for the bilinear coupling of equation 
(3.2) can easily be computed using the formula of equation ( 2.12f ). At one-loop, we 
find 



(Att) 2 (5^ =2M-Y*-Yj + 2Y U -Y}-M 
and the two-loop part of the /^-function is given by 

M 



(3.14) 



(^) 4 /i? 



9 V* vT y* yT 
^ 1 v 1 e 1 e ' 1 v 



9 V* y* yT 

* 1 V ' 1 V 1 V ' 1 V 



+ 



- 6 y; • yj Tr(y u • yj) - 2 y; • yj tt(y„ ■ yj) 

+ ^ 2 y;-yJ + 6 5 2 y;-yJ 
5 

2 Y y ■ Y e ■ Y e • Yj 2 Y v ■ Yj • Y v • yj" 

• 6 y, • yj Tr(y M • y„t) - 2 y, • yj Tr(y, • yj) 



+ - 5 2 y,-y 1 t + 6. 9 2 y,-yt 



52 



• M . 



(3.15) 



- 11 - 



In typical models of neutrino masses based on the see-saw mechanism, the effective 
neutrino mass operator of equation (|3.8| ) is obtained by integrating out the singlet 
superfields, which leads to the relation k = 2Yj 'M _1 Y^. The /^-function for M, 
together with the j3- function for Y v of equation ( 3.6d| ) , is therefore required to evolve 
the neutrino mass matrix from the GUT scale to the scale of the largest eigenvalue 
of M, if the singlets have a direct mass term. 



4. Discussion and Conclusions 



We have presented a general method to calculate two loop /3-functions for renor- 
malizable and non-renormalizable operators of the superpotential using superfield 
techniques. This method is very useful for model building since it provides a con- 
struction kit which allows to calculate the /^-functions in a given supersymmetric 
GUT model with little effort. We have applied this method to calculate the two- 
loop beta functions for the lowest-dimensional effective neutrino mass operator in the 
MSSM and for the the Yukawa couplings and the mass matrix in the MSSM extended 
by singlet chiral superfields. We have computed and specified the wavefunction re- 
normalization constants in the latter model, from which, using our method, the two 
loop RGE's for every, even higher dimensional, operator of the superpotential can 
directly be computed. A classification of several higher-dimensional operators for 
generating neutrino Majorana masses is e.g. given in [23]. Many of them can be gen- 
eralized for supersymmetric models. Their /^-functions can easily be obtained with 
our method. The two loop /3-function for the lowest-dimensional effective neutrino 
mass operator is required to increase the accuracy of many studies based on the one 
loop RGE [9—11], e.g. [24-32]. This accuracy may be needed for the neutrino sector 
since due to the absence of hadronic uncertainties, high precision measurements of 
the neutrino parameters may be achieved in future experiments. 
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